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Bulk viscosity of gauge theory plasma at strong coupling 
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We propose a lower bound on bulk viscosity of strongly coupled gauge theory plasmas. Using 
explicit example of the M = 2* gauge theory plasma we show that the bulk viscosity remains finite 
at a critical point with a divergent specific heat. We present an estimate for the bulk viscosity of 
QGP plasma at RHIC. 
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Recently, a holographic link between finite tempera- 
ture gauge theories and string theory black holes emerged 
as a viable theoretical tool to model properties of strongly 
coupled quark gluon plasma (QGP) produced at RHIC 
d 0> IE El- While the precise holographic dual to QCD 
is still missing, a progress in study of string theory black 
holes made it possible to compare the thermodynamics 
of strongly coupled QCD-like gauge theories [E ||| with 
lattice results |7|. The dual holographic approach has 
been successful to address dynamical properties of QGP 
such as the shear viscosity [8j and the parton jet quench- 
ing [EEli where few alternative techniques are available. 
Intriguingly, dual string theory studies reveal certain uni- 
versal features of gauge theory plasma dynamics. A no- 
table examples is the ratio of the shear viscosity^ to the 
entropy density s. It was shown in [D], HE EE 13 that 
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in any gauge theory plasma at infinite 't Hooft coupling, 
irrespectively of the dimensionality of the space, the mi- 
croscopic scales of the theory, and chemical potentials for 
the conserved quantities. The universality of the shear 
viscosity ratio ([T]) in strongly coupled gauge theories at 
finite temperature led Kovtun, Son and Starinets (KSS) 
to conjecture a shear viscosity bound (l5j 
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for all physical systems in Nature. Empirically, the KSS 
bound indeed appears to be satisfied by all common sub- 
stances [HI; moreover, it is correct at large (but finite) 
't Hooft coupling in J\f = 4 Yang-Mills theory plasma 

Em. 

We believe that it is such universal features of dual 
holographic models of gauge theories that might have 
some relevance to QCD. Thus, it is imperative to ask 
what are other generic properties of strongly coupled 
gauge theories. The question is complicated as neither 



the bulk viscosity [18[ nor the quenching of parton jets 
[l9| is universal for different gauge theory plasmas. 

It this Letter we propose a lower bound on bulk vis- 
cosity £ of strongly coupled gauge theories. Based on 
holographically dual computations, we conjecture that a 
bulk viscosity in a strongly coupled gauge theory plasma 
in p-space dimensions satisfies 
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where c s is the speed of sound. Notice that unlike the 
shear viscosity bound |(3J), our bound ([3]) is dynamical: 
as the temperature varies, generically both the speed 
of sound and the ration of bulk-to-shear viscosities will 
change. Our claim is that the bound ([3J is correct over 
all range of temperatures. 

In the following we present evidence in support of the 
bulk viscosity bound ©. First, we observe that the 
bound is saturated by the p + 1 space-time dimensional 
gauge theory plasma holographically dual to a stack of 
near-extremal flat Dp-branes [2(J, as well as in the hy- 
drodynamics of Little String Theory [IE EH- Second, 
we point out that the bound ^ remains saturated once 
above p-space dimensional gauge theory is compactified 
on a k < p space-dimensional torus [IE El ■ Third, we 
observe that the bound is satisfied (but in general not sat- 
urated) in certain 3 + 1 strongly coupled non-conformal 
plasma at high temperature [H, El . Finally, we present 
results for the bulk visco sity of the Af = 2* gauge 
theory plasma [E EE EE EE EE El over a wide range 
of temperatures, and for various mass deformation pa- 
rameters. We find that the bulk viscosity of the JV = 2* 
plasma satisfies the bound ([3]). As observed in [5j, the 
J\f = 2* plasma with zero fermion masses undergoes 
an interesting phase transition with vanishing speed of 
sound. A detailed analysis of the critical point 24} re- 
veals that at the transition point the specific heat di- 
verges as cy ~ |1 — T c /T| -1 / 2 . We find that despite the 
divergent specific heat the bulk viscosity at criticality re- 
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mains finite. We use results for the N = 2* gauge theory 
plasma to estimate the bulk viscosity of QGP at RHIC. 

Bulk viscosity of Dp-brane gauge theory plasma. M = 4 
Yang-Mills plasma at strong coupling is holographically 
dual to near-extremal stack of D3 branes. In this case 
conformal invariance of the theory implies that 
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Eq. ((¥]) was verified in supergravity approximation in [30| 
and beyond the supergravity approximation in [l7l |. No- 
tice that M = 4 plasma trivially satisfies the bound©. 

In [2(| the authors generalized computation of [30|to 
p + 1 space-time dimensional gauge theory plasma holo- 
graphically dual to near-extremal stack of Dp branes. 
They found the following dispersion relation for the 
sound waves 
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Hydrodynamics of a fluid with shear and bulk viscosities 
{?7,£} in p-space dimensions predicts the following sound 
wave dispersion 
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Using the universality of the shear viscosity (JTJ , one can 
verify that the bound ([3]) is saturated [2(| in the hydro- 
dynamics of the flat Dp branes. It is saturated as well in 
the hydrodynamics of Little String Theory [2(| HH • 

We point out now that the bound ([3]) is saturated as 
well for above strongly coupled gauge theory plasmas 
compactified on a /c-dimensional torus ( k < p ) [32| . 
Indeed, upon such a compactification the dispersion re- 
lation ([5]) will not change — much like an equation of 
state it is sensitive only to the local properties of the 
background geometry: 
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On the other hand, the hydrodynamics relation ([7]) is sen- 
sitive to the number of macroscopic ( infinitely extended 
) directions: 
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Again, using the universality of the shear viscosity |TJ) we 
find ( see also Eq. (5.2) of Ref. ) 
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It is precisely for the stated reason the bound is satu- 
rated in Sakai-Sugimoto model in the quenched approxi- 
mation [lH , even though 
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Bulk viscosity of non- conformal plasma at high tem- 
peratures. A much more nontrivial example is the bulk 
viscosity of non-conformal gauge theory plasma in four 
dimensions. The computation in the cascading gauge 
theory HH, |1] produced ^ 
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where A is the strong coupling scale of the cascading 
gauge theory. 

Likewise, for N —2* gauge theory plasma with bosonic 
and fermionic mass deformation parameters mj, <C T and 



1 



2 R T 

3 



16 



O 



3 s 



(13) 



where /?F « 8.001 [1 
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where (3j w 0.66666 

In all cases above we find that the viscosity bound © 
remains true — in general, it is no longer saturated. 

Bulk viscosity of M — 2* plasma. The strongest sup- 
port for the bulk viscosity bound ([3]) comes from study of 
the M = 2* bulk viscosity over the wide range of temper- 
atures. Such analysis is a direct extension of the frame- 
work presented in [l8| . The computations are extremely 
technical and will be detailed elsewhere (24|. Here, we 
report only the results of the analysis [35j . 

Fig. Q] represents the ratio * versus the speed of sound 
in M = 2* gauge theory plasma with rrif — 0. This model 
reaches a critical point with vanishing speed of sound at 
^ sa 2.32591 [5|. Although near the critical point the 

specific heat diverges as cy ~ |1 - T c /T| -1 / 2 [24] ( also 
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FIG. 1: Ratio of viscosities ^ versus the speed of sound in 
Af = 2* gauge theory plasma with zero fermionic mass defor- 
mation parameter ra/ = 0. The dashed line represents the 
bulk viscosity bound @. 





6.66 














6.64 














6.62 












-0.0005 




0.0005 


0.0010 


0.0015 


. 0.0020 



FIG. 2: Ratio of viscosities * in Af = 2* gauge theory plasma 
near the critical point. 



Fig. 8 of Q ), we find that the bulk viscosity remains 
finite, Fig. [2 and Fig. [3] . 

Fig. [4] represents the ratio ^ versus the speed of sound 
in M — 2* gauge theory plasma with "supersymmetric" 
mass deformation parameters — mj — m. We did 
not find any phase transition in this system up to tem- 
peratures T ss -j^. 

The dashed line in Fig.Q]and Fig.|4]represents the bulk 
viscosity bound ©• In both cases the bound is satisfied. 

Estimates for the viscosity of QGP at RHIC. It is 
tempting to use the Af = 2* strongly coupled gauge the- 
ory plasma results to estimate the bulk viscosity of QGP 
produced at RHIC. For c 2 s in the range 0.27 — 0.31, as in 
QCD at T = l.5T deconfinement [H, H3| we find 
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FIG. 3: Ratio of viscosities * in = 2* gauge theory plasma 
with zero fermionic mass deformation parameter m/ = 0. 
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FIG. 4: Ratio of viscosities a versus the speed of sound in 
M = 2* gauge theory plasma with "supersymmetric" mass 
deformation parameters mj, = TO/ = to. The dashed line 
represents the bulk viscosity bound ©. We computed the 
bulk viscosity up to m/T ~ 12. A single point represents 
extrapolation of the speed of sound and the viscosity ratio to 
T — » +0. 



Since RHIC produces QGP close to its criticality, we be- 
lieve that to/ = Af = 2* gauge theory model would re- 
flect physics more accurately. If so, it is important to re- 
analyze the hydrodynamics models of QGP with nonzero 
bulk viscosity in the range given by (|15p . 

In this Letter we presented some evidence in support 
of the bulk viscosity bound in strongly coupled gauge 
theory plasmas. It would be interesting to examine other 
holographic models and test the bound. As in [l2j], it 
would be interesting to study applicability of the bound 
in common substances realized in Nature. It appears 
that common liquids, like water, satisfy the bound [HI]. 
While the bound is generically satisfied in polyatomic 
gases 39], it is violated in monoatomic gases 40]. The 
bound also appears to be violated in high-temperature 
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QCD at weak coupling [4l[. In fact, experimental study 
of the bulk viscosity in argon at different densities [42| 
demonstrates that its ratio of bulk-to-shear viscosities 
violates/satisfies the bound at small/large densities. All 
this indicates the relevance of the bulk viscosity bound 
([3]) to strongly coupled systems only. 

We demonstrated that the bulk viscosity in the JV = 2* 
plasma with vanishing fermionic masses has a finite vis- 
cosity at the critical point with divergent specific heat. 
The corresponding critical exponent a = 0.5 ( cy ~ 
|1 — T c /T\~ a ) coincides with the mean-field universal 
value at the tricritical point [431 ] . Such a tricritical point 



is realized experimentally in solids [44j. It would be in- 
teresting to find a fluid with such a universal tricritical 
point and compare its bulk viscosity with that of the 
Af = 2* plasma at criticality. 
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